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Abstract 

We study the statistical mechanics of a one-dimensional log gas with general potential 
and arbitrary /3, the inverse of temperature, according to the method we introduced for 
two-dimensional Coulomb gases in |SS2j . Such ensembles correspond to random matrix 
models in some particular cases. The formal limit (3 = oo corresponds to "weighted Fekete 
sets" and is also treated. 

We introduce a one-dimensional version of the "renormalized energy" of [SSI] , mea- 
suring the total logarithmic interaction of an infinite set of points on the real line in a 
uniform neutralizing background. We show that this energy is minimized when the points 
are on a lattice. 

By a suitable splitting of the Hamiltonian we connect the full statistical mechanics 
problem to this renormalized energy W, and this allows us to obtain new results on the 
distribution of the points at the microscopic scale: in particular we show that configu- 
rations whose W is above a certain threshhold (which tends to min W as (3 — > oo) have 
exponentially small probability. This shows that the configurations have increasing order 
and crystallize as the temperature goes to zero. 

> 
OO 

vo 1 Introduction 

On 

In |SS2| we studied the statistical mechanics of a 2D classical Coulomb gas (or two-dimensional 
plasma) via the tool of the "renormalized energy" W introduced in |SS1| . a particular case 
of which is the Ginibre ensemble in random matrix theory. 

In this paper we are interested in doing the analogue in one dimension, i.e. first defining 
a "renormalized energy" for points on the real line, and applying this tool to the study of the 
classical log gases i.e. to probability laws of the form 



X 

(1.1) C/P£(X1,...,0 = J ?e -!™»(*l.".,*n) da . 1 ...dx r . 

where Zn is the associated partition function, i.e. a normalizing factor such that Pn is a 
probability, and 

n 

(1.2) w n (xi, ...,x n ) = - ^log|xj - Xj\ +n'Y]v(xi). 

Here the Xj's belong to R, /3 > is a parameter corresponding to (the inverse of) a tem- 
perature, and V is a relatively arbitrary potential, satisfying some growth conditions. For 
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a general presentation, we refer to the textbook |For| . Minimizers of w n are also called 
"weighted Fekete sets" and arise in interpolation, cf. |SaTo| . 

There is an abundant literature on the random matrix aspects of this problem (the con- 
nection was first pointed out by Wigner and Dyson |Wi[ |Dy| ) , which is the main motivation 
for studying log gases. Indeed, for the quadratic potential V{x) = x 2 /2, particular cases 
of (3 correspond to the most famous random matrix ensembles: for (3 = 1 the law Wn is 
the law of eigenvalues of matrices of the Gaussian Orthogonal Ensemble (GOE), while for 
(3 = 2 it corresponds to the Gaussian Unitary Ensemble (GUE), for general reference see 
[Fort IAGZ[ IMej . For V(x) still quadratic, general /3's have been shown to correspond to 
tri-diagonal random matrix ensembles, cf. |DE[ IABF] , This observation allowed Valko and 
Virag |VV] to derive the sine-/3 processes as the local spacing distributions of these ensembles. 
When (3 = 2 and V(x) is more general, the model corresponds (up to minor modification) 
to other determinantal processes called orthogonal polynomial ensembles (see e.g. |Koj for a 
review) . 

The study of P« via the random matrix aspect is generally based on explicit formulas for 
correlation functions and local statistics, obtained via orthogonal polynomials, as pioneered 
by Gaudin, Mehta, Dyson, cf. }Me[ lDj IDG| . We emphasize that for general (3 and V, 
which is the setting of the present work, there is no random matrix model associated and no 
explicit formulas available, and fewer results. Notable exceptions are the work of Ben Arous 
and Guionnet jBGj also valid for every (3 and the recent work of Bourgade, Erdos and Yau 
|BEYH lBEY2j . which are in the same setting as ours. We will discuss them more below. The 
global and local statistics (e.g. spacings) of eigenvalues of large random matrices are expected 
to be independent of the fact that the entries are specifically Gaussian, this is refered to as 
"universality", cf. |TV| IEPRSY| for latest results. Here we are rather interested in questions 
of universality with respect to the potential V, as in [B EY1| |BEY2| . 

The results here are counterparts of the results we obtained in |SS2j for xi,...,x n be- 
longing to M. 2 (this corresponds for V quadratic and fj = 2 to the Ginibre ensemble of non- 
hermitian Gaussian random matrices), or in other words the two-dimensional Coulomb gas. 
The study in |SS2| relied on relating the Hamiltonian w n to a Coulomb "renormalized en- 
ergy" W introduced in [SSI] in the context of Ginzburg-Landau vortices. This relied crucially 
on the fact that the logarithm is the Coulomb kernel in two dimensions, or in other words 
the fundamental solution to the Laplacian. When looking at the situation in one dimension, 
i.e. the present situation of the ID log-gas, the logarithmic kernel is no longer the Coulomb 
kernel, and it is not a priori clear that anything similar to the study in two dimensions can 
work. Note that the ID Coulomb gas, corresponding to where the logarithmic interaction 
is replaced by the ID Coulomb kernel \x\, has been studied, notably by Lenard [LelL ILe2] . 
Brascamp-Lieb [BL] , Aizenman-Martin [AM| . The situation there is rendered again more ac- 
cessible by the Coulomb nature of the interaction and its less singular character. In particular 
|BL] prove cristallization (i.e. that the points tend to arrange along a regular lattice) in the 
limit of a small temperature, we will get a similar result for the log-gas. 

The starting point of our study is that even though the logarithmic kernel is not Coulombic 
in dimension 1 , we can view the particles on the real line as embedded into the two-dimensional 
plane and interacting as Coulomb charges there. This provides a way of defining an analogue 
of the "renormalized energy" of [SSJ in the one-dimensional setting, still called W, which 



goes "via" the two-dimensional plane, see Section 1.1 below. 

Once this is accomplished, we connect in the same manner as [SS2] the Hamiltonian w n to 
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the renormalized energy W via a "splitting formula" (cf. Lemma 1.7 below), and we obtain 
the counterparts results to |SS2] : 

• a next-order expansion of the partition function in terms of n and f3, cf. Theorem |3j 
For V quadratic, exact formulas were already known via Selberg integrals. 

• the proof that ground states of w n , or "weighted Fekete sets", converge to minimizers 
of W, cf. Theorem |4j 

• the proof that the minimum of W is achieved by the one-dimensional regular lattice Z, 
called the "clock distribution" in the context of orthogonal polynomial ensembles [51] . 
This is in contrast with the dimension 2 where the identification of minimizers of W is 
still open (but conjectured to be "Abrikosov" triangular lattices.) 

• A large deviations type result which shows that events with high W become less and 
less likely as (3 — > oo, proving in particular the crystallization as the temperature tends 
to 0, in view of the result of the previous item. 

As far as we know, this is the first proof of crystallization for one-dimensional log gases. Our 
renormalized energy W, which serves to prove the crystallization also appears (like its two- 
dimensional version) to be a measurement of "order" of a configuration at the microscopic 
scale 1/n. What we show here is that there is more and more order (or rigidity) in the log gas, 
as the temperature gets small. Of course, as already mentioned, it is known that eigenvalues 
of random matrices, even of general Wigner matrices, should be regularly spaced. In fact the 
local statistics of the random process are completely known from |VVj . and [BEY1, BEY2 
showed that this could be extended to general Vs. Our results approach this question sort 
of orthogonally, by exhibiting a unique number which measures the average rigidity, without 
using explicit formulas for local statistics relying on determinantal forms. Note that in jBSel 
the second author and Borodin used W as a way of quantifying the order of random point 
processes, in particular those arising as local limits in random matrix theory. 

Our study here differs technically from the two-dimensional one in two ways: the first one 
is in the definition of W by embedding the problem into the plane, as already mentioned. The 
second one is more subtle: in both settings a crucial ingredient in the analysis is to reduce the 
evaluation of the interactions to an extensive quantity (instead of sums of pairwise Coulomb 
interactions); that quantity is essentially the L 2 norm of the "electric field" generated by the 
Coulomb charges. Test-configurations can be built and their energy evaluated by "copying 
and pasting" , provided a cut-off procedure is devised: it consists essentially in taking a given 
electric field and making it vanish on the boundary of a given box while not changing its 
energy too much. In physical terms, this corresponds to screening the field. The point 
is that screening is much easier in two dimensions than in one dimension, because in two 
dimensions there is more geometric flexibility to move charges around. The screening result 
of [SS2] does not apply to the 1-dimensional setting. We found that in fact, in dimension 1, 
not all configurations with finite energy can be effectively screened. However, we also found 
that generic "good" configurations can be, and this suffices for our purposes. The screening 
construction, which is different from the two-dimensional one, is one of the main difficulties 
here, and forms a large part of the paper. 

The rest of the introduction is organized as follows: in Section [l.l[ we present the embed- 
ding of the problem into two- dimensions and the definition of the one- dimensional renormal- 



ized energy W. In Section 1.2 we introduce notation and well-known facts on the equilibrium 
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measure (i.e. the minimizer of the mean- field limiting Hamiltonian), as well as the "splitting 



formula" that connects w n to W . Finally, in Section 1.3 we state our main results. 



1.1 Definition of the renormalized energy in 1 dimension 

As mentioned above, the renormalized energy between points in ID is obtained by "embed- 
ding" the real line in the plane and computing the renormalized energy in the plane, as defined 
in [SSlj . More specifically, we introduce the following definitions: 

M will denote the set of real numbers, but also the real line of the plane M 2 i.e. points of 
the form (x,0) S M 2 . For the sake of clarity, we will denote points in R by the letter x and 
points in the plane by z = (x,y). For a function x on K, we define its natural extension x 
to a function on IR 2 by xi x -> y) := x( x )- Ir wm denote the intervals [-R/2, R/2] in R and Ir 
the strips Ir x R in IR 2 . 5^ denotes the measure of length on M seen as embedded in IR 2 , that 
is 

ip5 R = / f(x, 0) dx 



for any smooth compactly supported test function tp in IR 2 . This measure can be multiplied 
by bounded functions on the real-line. For any measurable set U, \U\ denotes its Lebesgue 
measure, and if U is a finite set #U denotes its cardinal. 5 P denotes the Dirac mass at p. 

The renormalized energy W is an energy of an infinite configuration of points on the real 
line, screened by a uniform background charge, again on the real line. It will later be obtained 
as a limit as n — > oo of logarithmic interactions of systems of n particles/points. 

We first introduce the classes A m corresponding to infinite configurations on the real line 
with density m. 

Definition 1.1. Let m be a nonnegative number. Let E be a vector field in M 2 . We say E 
belongs to the admissible class A m if 

(1.3) div E = 2n{v - m5 R ), curl£ = in R 2 
where v has the form 

(1.4) v = S p for some discrete set AdC M 2 , 

pgA 

and 

(1.5) -j^- is bounded by a constant independent of R > 1. 

Definition 1.2. Let m be a nonnegative number. For any bounded function x an d an U 



E G A m satisfying (j 1 . 3|) we let 
(1.6) 



W(E, X ) = lim - / xl^l 2 + vrlog77 Vx(p) 

V^O J R 2\U p€A B(p, V ) ^ 



where x is the natural extension of X- 



4 



Here we thus look at any configuration of points in the plane through an "electric field" 
E that it generates in the whole plane M 2 . Observe that if E satisfies (1.3), there exists a 



function H such that E = — ViT, we will call this function the potential generated by the 
point charges. We also note that the distribution of charges on the real line is compensated 
by a background charge m5^ which is also concentrated on the real line. We will use the 
notation \R f° r positive cutoff functions over R satisfying, for some constant C independent 
of R, 

(1.7) |Vxh| < C, SuppOo?) C I R , X r{x) = 1 if \x\ < R/2 - 1. 
Definition 1.3. The renormalized energy W is defined, for E £ A m , by 

(1.8) W(E) = limsup W ( E ' XR \ 

R-^oo R 



where {xr}r>o satisfies (1.7). 



Returning to writing E = — V-ff, we see that while W in 2D can be viewed as a "renormal- 
ized" way of computing ||/f||^i( R 2), in ID it amounts rather to a renormalized computation 
of ||^i/2( K ) (where H s denote the fractional Sobolev spaces). In other words, because the 
logarithmic kernel is not Coulombic in one-dimension, the associated energy is non-local (and 
the associated operator is the fractional Laplacian A 1 / 2 ). Augmenting the dimension by 1 
allows to make it local and Coulombic again. This well-known extension idea is also now 
commonly used in the study of general fractional Laplacians [CSJ. 

As in the two dimensional case, we have the following properties: 



- The value of W does not depend on {xr}r as long as it satisfies (1.7). 

- W is insensitive to compact perturbations of the configuration. 

- Scaling: it is easy to check that if E belongs to A m then E' := —E(-/m) belongs to A\ 
and 

(1.9) W{E) = m(W(E') -Trlogm) . 

- If E £ A m then in the neighborhood of p G A we have div E = 2ir(5 p — md^.), curli? = 0, 
thus we have near p the decomposition E(x) = —V log \x — p\ + f(x) where / is smooth, 



and it easily follows that the limit (1.6) exists. It also follows that E belongs to Lf oc for 
any q < 2. 

An extra fact proven below is valid only in the one-dimensional case (and not the two- 
dimensional one): 



Lemma 1.4. Let E £ A m be such that W(E) < +oo. Then any other E' satisfying (1.3) 



(1.4) with the same v and W(E') < +oo, is such that E' = E. In other words, W only 



depends on the points. 

By simple considerations similar to [SS2, Section 1.2] this makes W a measurable function 
of the bounded Radon measure v. 
This implies in particular 
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Corollary 1.5. Under the same assumptions, if S(x, y) = (x, —y) then E o S = S o E. 



S o E o S satisfies (1.3) with the same v as E, and 



Indeed, it is easy to check that E' 
obviously W(E') < +00, hence E' = E. 

The following lemma is proven in [BScJ, and shows that there is an explicit formula for 
W in terms of the points when the configuration is assumed to have some periodicity. Here 
we can reduce to m = 1 by scaling, as seen above. 

Lemma 1.6. In the case m = 1 and when the set of points A is periodic with respect to 
some lattice NZ, then it can be viewed as a set of N points a\, . . . , a/y over the torus Tjy := 



/(iVZ). In this case, by Lemma 1.4 there exists a unique E satisfying (1.3) and for which 



W(E) < +00. It is periodic and equal to Er a .\ = VH, where H is the solution on TV to 
—AH = 2n(^2i 5 ai — Sjh), and we have the explicit formula: 



(1.10) 



W(E {ai} ) = -^J2 l °' 



2 sin 



7T (Oj 



N 



IT fog 



2tt 
77' 



As in the two-dimensional case, we can prove that min v 4 m is achieved, but contrarily to 
the two-dimensional case, the value of the minimum can be explicitly computed: we will prove 
here the following 



Theorem 1. min^ 
i.e. A = iZ. 



W 



-Tim log(27rm) and this minimum is achieved by the perfect lattice 



We recall that in dimension 2, it was conjectured in jSSlj but not proven, that the min- 
imum value is achieved at the triangular lattice with angles 60° (which is shown to achieve 
the minimum among all lattices), also called the Abrikosov lattice in the context of super- 
conductivity. 

The proof of Theorem [T] relies on showing that a minimizer can be approximated by 
configurations which are periodic with period — > 00 (this result itself relies on the screening 
construction we mentioned above), and then using a convexity argument to find the minimizer 



among periodic configurations with a fixed period via (1.10). 



1.2 The equilibrium measure and the splitting formula 

The Hamiltonian ( |1.2[ ) is written in the mean-field scaling. The limiting "mean-field" limiting 
energy and its minimizer are in fact well-known, and they tell us the optimal macroscopic 
distribution of points — what we investigate here is the microscopic distribution of points, 
as a correction to this optimal macroscopic behavior. More precisely, the mean- field limit 
energy is 

(1.11) F{n) = j — log I a; — y\ d/j,(x) dfi(y) + / V(x)d/j(x), 

JRxR ' ' Jr 

and its unique minimizer among probability measures is called the equilibrium measure in the 
language of potential theory (cf. [SaTo| ) . We will denote it //o- It is not hard to prove that the 
"spectral measure" (so-called in the context of random matrices) v n = - X^iLi °~Xi converge to 
/io, the minimizer of J 7 , itself also called Voiculescu's noncommutative entropy in the context 
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of random matrices, cf. e.g. |AGZ] and references therein. The sense of convergence usually 
proven is 



(1.12) FP (V/ E C b (C,R),J fdv n ^ J fdfio 

For example, for the case of the GUE i.e. when V(x) = \x\ 2 and (3 = 1, the correspond 
distribution fj,Q is simply Wigner's "semi-circle law" p(x) = ^ \/4 — x 2 l|a,| <2 , cf. [Wil IMe| . A 
stronger result was proven in |BG| for all (5 (cf. [AGZj for the case of general V): it estimates 
the large deviations from this convergence and shows that T is the appropriate rate function. 
The result can be written: 

Theorem 2 (Ben Arous - Guionnet [BGj ). Let f3 > 0, and denote by Fn the image of the law 
( 1.1 ) by the map {x\, . . . , x n ) i->- v n , where v n = j i Y^t=\ $xi ■ Then for any subset A of the set 



of probability measures on R (endowed with the topology of weak convergence), we have 
- inf T{p) < liminf -= logP^(^) < limsup \ logP^(A) < - inf T(n), 

where T = T — min T . 

Our result Theorem [5] below will be a sort of next-order extension of this result. Note that 
the Central Limit Theorem for (macroscopic) fluctuations from the limiting law was proved 
by Johansson [Jo] . 

Let us now state a few facts that we will need about the equilibrium measure (io, for which 
we refer to [SaToj: is characterized by the fact that there exists a constant c (depending 
on V) such that 
(1.13) 

V V 
U^ + — = c quasi-every where in the support of [io, and + — > c quasi-every where 

where for any /x, is the potential generated by \x, defined by 



(1.14) U fl (x) = - / log \x - y\ dfi(y). 

(Note that looking in Fourier, this corresponds to = — A -1 / 2 /^). Here and in all the paper, 
we denote by A -1 the operator of convolution by ^ log | • | in M 2 . It is such that Ao A -1 = Id, 
where A is the usual Laplacian. We may view /xo as a measure over M 2 and extend into 
a function of R 2 by U»(z) = - f R log\z - (y,0)\ dfi(y). Then If 1 = -A _ V in R 2 . 
We also define 

V 

( = TJM + — - c 



where c is the constant in (1.13). From the above we know that ( > in R and £ = in 
E := Supp(po)- We will make the assumption that /xo has a density mo with respect to the 
Lebesgue measure. 

It's now time to state our assumptions on V: we assume in the sequel 

V(x) 

(1.15) V is lower semicontinuous and lim — log|sc| = +oo. 

\x\— >+oo 2 
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(1.16) 



S is a finite union of closed intervals Si, ... , Ej^. 



(1.17) There exist 7, m > such that j\/ dist (x, M. \ S) < mo(x) < m for all x 6 



(1.18) 



m G C°'5(M). 



(1.19) 



There exists /3i > such that 



-ftCV^-loglsl)^ 



< +00. 



-1,1 



The assumption ( 1.15 ) ensures (see [SaTo| ) that ( 1.11 ) has a minimizer, and that its support S 
is compact. Assumptions ( 1.16 )-( 1.18) are needed for the construction in Section|4j They are 
certainly not optimal but are meant to include at least the model case of po = P, Wigner's 
semi-circle law. Assumption (1.19) is a supplementary assumption on the growth of V at 



infinity, needed for the case with temperature. It only requires a very mild growth of V/2 
log|x|, i.e. slightly more than (1.15). 



We are interested in examining the "next order" behavior, or that of fluctuations around 
the limiting law fj®. While J- and po are found through the large deviations at speed n 2 , as seen 
in Theorem [2| we look into the speed n and show, not a large deviations principle, but that a 
suitably defined average J W dP of the "renormalized energy" W defined above (1.8) acts as 
a sort of rate function, and that a "threshhold phenomenon" holds: JWdP < miniy + Cp 
except with exponentially small probability. The crucial fact is that Cp — > as (5 — > 00, which 
shows the crystallization. The average of W corresponds to the average of W computed over 
all blown-up configurations of the points (with respect to all possible blow-up centers in S) 
at the scale n. The measure P is like a Young measure on all the possible limits of blow-ups 
of E. 

The starting point of our analysis is to establish, through a simple but crucial exact 
"splitting formula", a link between the law and the energy w n , and W: 

Lemma 1.7 (Splitting formula). For any x\, . . . , x n € M the following holds 

1 n 
(1.20) w n (x u ...,x n ) = n 2 T(fi ) ~ n log n + -W(VH' n , 1 R2 ) + 2nV(( 

7T ' — ' 



X 1 



1=1 



where H' n is the solution in M 2 to 



(1.21) 



H 'n( x ') = -27rA 1 I ^ S x > ; - m (x'/n)5M 



,i=l 



with x\ = nxi (viewed as points in M?), and W(-, •) is as in (1.6) 
We may then define 



(1.22) 



F n {y) 



i(iW(Vir; j l R9 ) + 2n/ R C^) if v is of the form E? =1 ^ 
+00 otherwise 



and also 



(1.23) 



F n (y) = F n (y)-2 / $dv<F n {v 
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and we thus have the following rewriting of w n : 

(1.24) w n (xi, ...,x n ) = n 2 J r ( / u ) - nlogn + nF n (v). 

This allows to separate orders in the limit n — > oo since one of the main outputs of our analysis 
is that F n {v) is of order 1. Moreover, F n (y) contains information about the microscopic 
structure of the configurations, its limit as n — > oo will be the average of W over all blown-up 
points, f W dP mentioned above. 

1.3 Main results 

The first main result is the announced one on the next order asymptotic expansion of the 
partition function, which becomes sharp as (5 — > oo. 



Theorem 3. Let V satisfy assumptions (1.15) — (1.19). There exist functions fx, fi depending 
only on V , such that for any (3q > and any (3 > (3q, and for n larger than some uq depending 
on (3q, we have 

(1.25) npfttf) < logZf - (-^n 2 F( m ) + ^nlogn) < n/?/ 2 (/3), 
with /i,/2 bounded in [/3o,+oo) and 

(1.26) Jim h{fi)= Urn fM = ~ 
where 

(1.27) a= ^o(x) log(27r^o(^)) dx. 



Remark 1.8. In fact we prove that the statement holds with /2(/3) = § + § for any C > 
log|E|. 



The exact value for Zn is only known for arbitrary /3's for V(x) = x 2 /2 via Selberg 
integrals cf. \M.e\ IAGZ] . For some more general potentials there are formulas for (3 = 2 (e.g. 
[EML, PAS, GMSj). For general potential and /3, this improves on the known results, which 
only gave the expansion log Zn ~ j3n 2 T(/Jo). 

Our next result is the one we mentioned about minimizers of w n , or weighted Fekete sets. 



It identifies the T-limit (in the sense of T-convergence, cf. [BrJ) of {F n } n , defined in (1.22) or 



(1.24). This allows a description of the minimizers at the microscopic level. Below we abuse 
notation by writing v n = X^=i ^» wnen it should be v n = Y27=i „ • For such a v n we let 
v' n = Ya=i fix', be the measure in blown-up coordinates x' = nx, and E Vn = —VH' n , where 



H' n is defined by (1.21). (To avoid confusion, we emphasize here that v n lives at the original 
scale while E Vn lives at the blown-up scale.) We set mo'(x') = mo(x'/n), where mo is the 
density of /io- We also let 



(1-28) P Vn = I 



fi(x,E Vn (nx+-)) dx, 



i.e. the push-forward of the normalized Lebesgue measure on X by x t— > (x, E Un (nx + •)). It 
is a probability measure on S x L^ oc (Mp ,M?) (couples of (blow-up centers, blown-up electric 
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field around this center)). In taking the T-limit of F n , the limiting object is more complex, 
it is the limit P of P„ n , i.e. a Young measure akin to the Young measures on micropatterns 
introduced in [AM]. Note that in the rest of the paper, the probability P has nothing to 
do with Fn, and depends on the realizations of the configurations of points. We state the 
generalized invariance property satisfied by these limiting object as a definition. 

Definition 1.9. We say a probability measure P on £ x L^ oc (M. 2 , M?) is T\( x yinvariant if P 
is invariant by (x, E) i— > (x, E(X(x) + •)), for any X(x) of class C 1 from £ to K. 

Theorem 4 (Microscopic behavior of weighted Fekete sets). Let the potential V satisfy as- 



sumptions (1.15)-(i-i8). Fix 1 < p < 2 and let X = £ x L? C 



2 )- 



A. Lower bound. Let v n = Y17=i a sequence such that F n (y n ) < C. Then P v 



defined by (1.28) is a probability measure on X and 



1. Any subsequence of {P Un }n has a convergent subsequence converging to a probability 
measure on X as n —> oo. We denote by P such a limit. 

2. The first marginal of P is the normalized Lebesgue measure on S. P is T\( x y invariant. 

3. For P almost every (x, E) we have E E A mQ ^ x y 
4- It holds that 



(1.29) 



liminf^K) > ^ / W(E)dP(x,E) > a, 

ri— >oo 7T 



where a is as in (1.27). 



B. Upper bound construction. Conversely, assume P is a T^m -invariant probability 
measure on X whose first marginal is j^r^|s an d such that for P-almost every (x,E) we 
have E E A mo ( x y Then there exists a sequence {v n = 2iLi$B»}n °f empirical measures on 
S and a sequence {E n } n in L^ C (M 2 ,M 2 ) such that div E n = 2ir(i/ n — and such that 

defining P n as in (1.28), with E n replacing E Vn , we have P n — > P as n 

(1.30) limsupF n (i/„) <— [ W{E)dP{x,E). 



co and 



C. Consequences for minimizers. If {x\, . . . ,x n ) minimizes w n for every n and v n 



Y17=i^^i> th> en the limit P of P Vn as defined in (1.28) satisfies the following. 

1. For P-almost every (x, E), the electric field E minimizes W over A mo ( x y 

2. We have 

lim F n {v n ) = lim T n {v n ) = M f W{E) dP(x, E) = a, lim V C(^) = 0. 

n— >oo n— ¥oo tx J n— >oo 

8=1 
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Note that part B of the theorem is only a partial converse to part A because the constructed 
E n need not be curl-free, hence in general E n ^ E Un . 

It can be expected that £ (which is positive exactly in the complement of S) controls 
the distance to £ (to some power). This can probably be proven under some strict growth 
assumption on V and explicit formulas for the Poisson kernel of the half-Laplacian. We do 
not pursue this question here. 

The next result is the announced next order result on deviations, to be compared to 
Theorem [2j 

We let A n be a subset of M n . We identify points in W 1 with measures v n of the form 
EiLi<W We also embed W 1 into the set of probabilities on X = £ x Lf oc (M 2 ,IR 2 ) in the 
following way: For any n and x = {x\, . . . , x n ) £ W 1 we let i n (x) = P Un , where u n = Y^i=\ 
and P Vn is as in ( 1.28| ), so that i n (x) is an element of V(X), the set of probability measures 



on X = £ x L^ oc (R 2 , IR 2 ), which we consider endowed with the topology of weak convergence. 



Theorem 5. Let V satisfy ( 1.15 ) — ( 1.19 ) . For any (3q > and any (3 > /3q the following holds. 



A. For any n > 0, let A n C M. n . Denote 



(1.31) A 00 =f]\J i m (A m ). 



n>0 m>n 



Then for any rj > there is C v > depending on V and rj only such that a being as in ( 1.27), 
logP^AO ^ /3/|S| . . [rrn^iv, „s Cv 



(1.32) limsu P logP " (An) <-^f^i inf [ W(E)dP(x,E) - a-r, 

n->oo n 2 V 7r PeA^ J 



Conversely, let A C V{X) be a set of T\( x y invariant probability measures on X and let A 
be the interior of A for the topology of weak convergence. Then for any n > 0, there exists a 
sequence of subsets A n C T, n such that 



(1.33) - - ( ^ inf f W(E)dP(x,E) -a + n+ ^] < lim inf 



logK(A n ) 



n 



and such that for any sequence {v n = X^ILi $xi\n such that (x%, . . . ,x n ) £ A n for every n 
there exists a sequence of fields E n € Lj 3 oc (IR 2 , M 2 ) such that div E n = 2n(v' n — and 



such that — defining P n as in (1.28) with E n replacing E v — we have 



(1.34) limP n £A. 

n 

B. For each integer n, and given j3 > (3q letQn denote the push-forward ofPn by i n defined 
above. Then Qn is tight in V(V(X)) and converges, up to a subsequence, to a probability 
measure Qr G V(V(X)). 



Note that if P n was Py n , then ( |1.34[ ) would be equivalent to saying that n n U m > n i (A m ) C A. 
The difference between P Vn and P n is that the latter is generated by a field E n which is not 
necessarily curl free (we believe this gap cannot be bridged). 

Part B of the result shows the existence of a limiting "electric field" process, hence equiv- 
alently, via projecting by (x,E) >->■ ^div E + mo{x)5K, of a limiting point process. 
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As we said, Part A of this theorem means that events such that the average of W is larger 
than a + r\ + C„//3 have exponentially decaying probability as n — > oo, with r\ > arbitrary. 
Note here the important fact that 



(1.35) min ^ J W(E) dP{x, E) = a 



where the min is taken over P probabilities satisfying P-a.e. (x,E) satisfies E G ■A rjlQ i x )i a 
fact which simply follows from (1.9) and Theorem [TJ 



This means that configurations should have more and more order (evaluated via W) as j3 
gets large, and have to crystallize to minimizers of W when /3 — > oo (the term crystallization 
can be legitimately used in view of Theorem [TJ . Note that we do not claim that this is a true 
large deviations result, since the corresponding lower bound is missing. In fact it is likely 
that the large deviations rate function is not truly W but contains an extra entropy term. To 
demonstrate its existence is an open question. 

Finally, let us mention that our method yields estimates on the probability of some rare 
events, typically the probability that the number of points in a microscopic interval deviates 
from the number given by ^o- We present them below, even though stronger results are 
obtained in [BEY1, BEY2] . The results below follow easily from the estimate (provided by 
Theorem 3| that F n < C except on a set of small probability. 



Theorem 6. Let V satisfy assumptions (1.15)-(1.18). There exists a universal constant 



Rq > and c, C > depending only on V such that: For any /3q > 0, any (3 > /3q, any n 
large enough depending on /3q, for any xi, . . . ,x n £ R, any R > Rq, any interval I C K of 
length R/n, and any rj > 0, letting v n = Y11=i we have the following: 

(1.36) logP^(K(I) -nn (I)\ >rjR) < -c^mm(n 2 ,7] 3 )R 2 + C/3(R + n) + Cn, 



(1.37) logP^ ((1 + R 2 /n 2 ) 1 * ^\vn-n^ \\ w -i, q{I) >n^j<-cnPn 2 + Cn(P+l) 

where W~ 1,q (I) is the dual of the Sobolev space Wq ,(1 (I), with X/q + 1/q' = 1, in particular 
W^ 1 ' 1 is the dual of Lips chitz functions; and 

(1.38) logWy Cdu n >r^j <-± n f3ri + Cn((3 + l). 
Note that in these results R can be taken to depend on n. 



(1.36) tells us that the density of eigenvalues is correctly approximated by the limiting 



law hq at all small scales bigger than Cn 1 / 2 for some C. However this in fact should hold at 



all scales with R > 1, cf. [ESYl IBEYlj IBEY2] , Q serves to control the probability that 

points are outside E (since {Q > 0} = S c ). 

The rest of the paper is organized as follows. In Section [2] we prove the preliminary results 
Lemma 1.4 and 1.7 as well as Theorem [TJ In Section [3] we turn to the core of the proof, with 



the statement of the main screening result and construction. Sections [4] and [5] contains the 
proofs. Finally Section [6] gathers the proofs of all the remaining results which are adapted 
from [SS2] . 

Acknowledgements: We are grateful to Alexei Borodin, Gerard Ben Arous, Amir Dembo, 
Percy Deift, Nicolas Fournier, Alice Guionnet and Ofer Zeitouni for their interest and helpful 
discussions. E. S. was supported by the Institut Universitaire de France and S.S. by a EURYI 
award. 
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2 The first results 



2.1 Proof of Lemma 11.41 



Assume E and E 1 belong to A m and satisfy (1.3) with the same v. Then / = E — E' is 
divergence- free and curl- free, hence can be seen, identifying M? and C, as an entire holo- 
morphic function Ylri=o a n zn - If we assume that W(E) and W(E') are finite, then it fol- 
lows from jSeTi] . Corollary 1.2 that for p < 2, the growth of the L 1 norms of E and E' 
is no worse than R 3 ^ 2 \/log R hence there exists C > such that for any R > 2 we have 
L^Bjf ) < CR?I 2 "v/l°g -R- But by Cauchy's formula we have, for any i? > and t G [i?, i?+l] 



It follows with the above that \a n \ < CR?l 2 yJ\og RR~ n ~ l which implies, letting R — > oo that 
a n = for any n > 1, thus / is a constant. This constant must then be zero since both E 
and ii/ are square integrable on the infinite strips [a, b] x [1, +oo]. 



2.2 Proof of the splitting formula (Lemma 1.7) 

Let v n = "522=1 Firet.) letting A denote the diagonal of M x 1, we may rewrite w n as 

w n (xi,...,x n ) = / - log \x - y\ dv n (x) dv n (y) + n / V(x) dv n (x). 
Splitting v n as ti^iq + v n — nno and using the fact that /io x Mo (A) = 0, we obtain 
w(x\, ... ,x n ) = n 2 J"(^ ) + 2n / U^°(x) d(y n - nfj, )(x) + n V{x) d{v n - n/j, )(x) 



+ / - log \x - y\ d(y n - n/jL )(x) d(v n - nfi Q )(y). 



Since + ^ = c + ( and since v n and n\iQ have same mass n, we have 
2n J U^°(x) d(u n - nfi )(x) + n J V(x) d(y n - n^)(x) = 2n J £ d(y n - /i ) = 2n J(dv n , 

using the fact that £ = on the support of //o- 
In addition, we have that 

(2.1) / - log \x - y\ d(y n - nfi )(x) d(u n - n/j, )(y) = -W(VH n , 1 R2 ), 

J(RxR)\A ^ 

where we define H n = — 2nA~ 1 (522=1 ^ — n ^o)- Indeed, the integral might as well be 
written as over M. 2 x M 2 \A with the diagonal in M. 2 xR 2 ; and then the identity is proven in 
[SS2, Section 2]. Combining all the above we find 

(2.2) w(xi, ...,x n )= n 2 F{^) + 2n j Qdv n + ^-W(VH n , 1 R2 ). 



13 



But, changing variables, we have 



\U? =1 B(xi, V ) 



R 2 \UY =l B(x' v n V ) 



iv<i 2 , 



and by adding irnlogrj on both sides and letting 77 — > we deduce that V7(V-ff n ,lj 
W(— VH' n , 1 K 2) — 7rnlogn. Together with (2.2) this proves the lemma. 

2.3 Minimization of W: proof of Theorem [I] 



By scaling (cf. (1.9)), we reduce to m = 1. The result relies on the fact that there exists a 
minimizing sequence for min_4 1 W consisting of periodic vector-fields: 

Proposition 2.1. For any large enough integer R, there exists a sequence {-ErIkgn in <A\ 
such that each Er is 2R-periodic (with respect to the x variable) and 

limsup W(Eji) < minW. 

The proof uses the main results of Section [3] and is postponed to the end of Section |6j 
The following proposition could be proven as in |SSlj . however we omit the proof here. 



Proposition 2.2. W : Lf oc (M 2 ,R 2 ) 
and is finite. 



U {+00} is a Borel function. inf_4 x W is achieved 



The result of Theorem [T] will follow from Proposition 2.1 combined with the following 

Proposition 2.3 (Minimization in the periodic case). Let a\, ■ ■ ■ , be any points in [0, N] 
and E{ a .y be the corresponding periodic vector field, as in Lemma 1.6. Then 

W(E {a .y) >W(E Z ) = -7Tlog27T 

where Ex is the electric field associated to the perfect lattice X. 



Proof. W(E%) is immediately computed via (1.10), taking N = 1. 

Let us now consider arbitrary points ai, . . . , ojy in [0, N], and assume a\ < ■ ■ ■ < ajy. Let 



us also denote = ai+i — ai, with the convention on+i = ai+N. We have 2~2^=i M M = N. 
Similarly, let u v ^ = ai +p — ai, with the convention a^v+z = ai + N. We have Yli=i u p,i = P-N- 
By periodicity of sin, we may view the points a^ as living on the circle K/ (NX). When adding 
the terms in ai — aj in the sum of ( 1.10 ), we can split it according to the difference p = j — i 
but modulo N. This way, there remains 
(2.3) 

> lncr "sin 

N 



W(E {ai} ) 



E 1 ^ 



N 



-7T log 



2tt 
~N 



2vr 

"N 



[N/2] N 

EE'-* 

p=i i=i 



2 sin ■ 



7TU 



p. 1 



T 



-7rlog 



2vr 

iV' 



where [•] denotes the integer part. But the function log |2sinx| is stricly concave on [0, ir]. It 
follows that 



1 

N 



N 



2 sin 



N 



< log 



N 



2 sin 



i=i 



log 



2 sin 



pir 
N~ 
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with equality if and only if all the u p i are equal. Inserting into (2.3) we obtain 

[N/2] 



(2.4) 



W{E {ai} ) > -27T lo § 
p=i 



2 sin 



pir 
~N 



7T log 



2vr 
iV' 



On the other hand, if we take for the o^'s the points of the lattice Z viewing them as N- 



periodic, we have u Pj i = p for all p, i, so if we compute W{E%) using (2.3), we find 

[N/2] N _ [N/2] 



W(E I 



2ir 

'~N 



EE lo f 

p=l i=l 



2 sin 



up 
~N 



7rlog 



2tt 
iV 



-2vr ^ log 
P =i 



2 sin 



7rp 
~N 



IT log 



2vr 
iV' 



This is the right-hand side of (2.4), so (2.4) proves that W(E{ a .y) > W{E-£) with equality if 
and only if all the equal, which one can easily check implies that {a^} = Z. □ 

Combining with Corollary |2.1[ this proves Theorem [T] 



3 The main screening result and construction 

This section contains the main "upper bound result" relying on explicit constructions, and 
requiring as a main ingredient the "screening result" mentioned in the introduction. We state 
the results, whose proofs will occupy the next section. 



Proposition 3.1. Let P be a T\i x y invariant probability measure (see Definition 1.9) on 



X = E x Ll 



loc V-"*- > - 



) with first marginal dxv^/\E\ and such that for P almost every (x, E) 
we have E £ ■A rm f x y Then, for any rj > 0, there exists 5 > and for any n a subset A n C W 1 
such that \A n \ > n\(5/n) n and for every sequence {v n = Y^i=i^m\n (z/ij • • • >3/n) £ A ri 
the following holds. 

i) We have the upper bound 

(3.1) limsup- (w n (y y n ) -n 2 J^(fi ) + nlogn) < ^ ( W{E) dP{x,E) + n. 

n— >oo n it J 

ii) There exists {E n } n in L^ oc (M 2 ,M 2 ) such that div E n = 2Tr(v' n — mo'S^_) and such that 
the image P n of dx\j^/\H\ by the map x \-t (x,E n (nx + •)) is such that 

(3.2) lim sup dist (P n , P) < n, 

n— >oo 

where dist is a distance which metrizes the topology of weak convergence on V(X). 

Applying the above proposition with rj = 1/k we get a subset A n ^ in which we choose 
any n-tuple (yi,k)i<i<n- This yields in turn a family {P n ,k} of probability measures on X. A 
standard diagonal extraction argument then yields 

Corollary 3.2 (Theorem[4j Part B). Under the same assumptions as Proposition 3.1, there 
exists a sequence {u n = Yl?=i $xi}n and a sequence {E n } n in L^ oc (IR 2 , M 2 ) such that div E n = 
27r(i/ — mo'dft) and 



(3.3) limsup- (w n (xi, ... ,x n ) 

n— >oo Tl 



n 



! T{fM ) + nlogn) < ^ j W{E)dP(x,E). 



Moreover, denoting P n the image of dx|s/|E| by the map x \— > (x, E n (nx + ■)), we have 
P n — > P as n — > +oo. 
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(3.4) 

where 
(3.5) 



We may cancel out all leading order terms and rewrite the probability law (1.1) as 

1 



n P 



X\ , . . . , x r 



K 13 = Z^e^ ( n2 - F (w))- nl °g ?1 ). 



l 2 F ^ dxi... dx r . 



A consequence of Proposition 3.1 is, recalling (1.27) and (1.35): 



Corollary 3.3 (Lower bound part of Theorem [3]). For any rj > there exists C v > such 
that for any (3 > we have 



(3.6) 



liminf — - — — > — — (0 + 77) — C v . 

n-H-00 n 2 



Proof of the corollary. It is exactly the same as in |SS2| but just letting 
(3.7) cr m E(y) := mE{my). 



□ 



We now state our main screening result, on which the proof of Proposition 3.1 is based, 
and which is the main difference with the two-dimensional situation : 



Proposition 3.4. Let Ir = [-R/2,R/2], let xr satisfy ( |L7|) . 

Assume G C A.\ is such that there exists C > such that for any E € G and writing 
we have 



v = ^-div E + 



(3.8) 



Vi? > 



1. sew 



R 



(3.9) 



lim W ^ XR) =W(E)<C, 



(3.10) lim lim f \E 



I f 

J Ir J \y\>ya 



2 



DO— H-00 R— >+oo 

and such that moreover all the convergences are uniform w.r.t. E G G. 

Then there exists for every < e < 1, there exists Rq > such that if R > Rq with 
R E N, then for every E E G there exists a vector field Er E L^ oc (I^,IR 2 ) such that the 
following holds: 

i) Er ■ v = on OIr, where v denotes the outer unit normal. 

ii) There is a discrete subset A C Ir such that 

div E R = 2ir ^ S p - <5 R in I R . 
Ipsa / 
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Hi) E R (z) = E(z) for x G [-12/2 + eR, R/2 - eR) . 



iv) 



(3.11) 



W(E R ,l lR ) 
R 



< W(E) + Ce. 



The assumption (3.10) is a supplementary assumption which is not satisfied for all vector- 



fields in Am-, an d which allows to perform the screening. We have not been able to show that 
screening is always possible without it. However we will see in Section [4] and in particular in 
Lemma |5.4| that this assumption is satisfied "generically" i.e. for a large set of vector-fields 
in the support of any invariant probability measure, and this will suffice for our purposes. 



4 Proofs of Propositions 3.4 and 3.1 



4.1 Preliminaries: a mass displacement result 

In this subsection, we state the analogue in ID of Proposition 4.9 in |SSlj . a result we will 
need later. The point of this result is to say that even though the energy density associated 
to W: \ \E\ 2 + 7r log r\ ]T) S p , is not positive or even bounded below, it effectively behaves as if 
it were, because it can be replaced by an energy-density g which is pointwise bounded below, 
at the expense of a negligible error. The density g is obtained by displacing the negative part 
of the energy- density into the positive part. The proof is no different than in [SSI] once the 
one-dimensional setting has been embedded into the two-dimensional one as stated. 

Proposition 4.1. Assume (f,E) are such that v = 27r^ pgA 5 p for some finite subset A of 
R, div E = 2tt(u — a(x)5^), for some a G L°°(M.), and cwclE = 0. Then, given < p < po, 
where po > is universal, there exists a measure density g in M? such that 

i) There exists a family of disjoint closed balls B p centered on the real line, covering 
Supp(u), such that the sum of the radii of the balls in B p intersected with any segment 
of R of length 1 is bounded by p and such that 



(4.1) 



9>-C(\\a\ 



L°° 



+ D + 



\E\ 2 1 



in 



where C depends only on p. 



ii) 



1 



\E\ in the complement of 



■1,1]. 



Hi) For any function x compactly supported in M we have, letting x( x ,y) = x( x )> 



(4.2) 



W(E,x)- / xdg 



<C#(togJV + ||a|| M )||Vx||oo, 



where N = 6 A | B(p, A) Pi Supp(X7x) ^ 0} and A depends only on p. 
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Proposition 4.9 of |SS1| of which the above Proposition is a restatement, was stated for 
a fixed universal po, but we may use instead in its proof any < p < po, which makes 
the constant C above depend on p. Another fact which is true from the proof of \SS1\ 
Proposition 4.9] but not stated in the proposition itself is that in fact g = ^\E\ 2 outside 
D p B(p,r) for some constant r > depending only on p, and if p is taken small enough, then 
we may take r = 1, which yields item ii) of Proposition 4.1 



The next lemma shows that a control on W implies a corresponding control on J g and 
of f \E\ 2 away from the real axis, growing only like R. 



Lemma 4.2. Assume that G C A\ is such that (3.8) and (3.9) hold uniformly w.r.t. E 6 G, 
for a certain constant Cq. Then for any E € G, writing v 
enough depending on G, we have 



2tt 



div E + 5m., for every R large 



(4.3) 



Wr)-R\ <C!ii 3 / 4 log R, 



(4.4) 



I R x{\y\>l} 



\E\ 2 < CR(W(E) + 1) 



and denoting by g the result of applying Proposition 4-1 to E for some fixed value p < 1/8, 
we have 



(4.5) 



W(E, X r) -Cii? 3 / 4 log 2 R< I dg < W(E, XR+i) + CiR 3 ^ log 2 R, 



where xr satisfies (1.7), C\ depends only on G and C is a universal constant. 

Proof. We denote by C\ any constant depending only on G, and by C any universal constant. 
From Q we have for any E G G that v{I R ) < CiR and W{E, X r) < C X R if R is 

large enough depending on G. Thus, applying (4.2) we have 

XRdg 



/ 



< CiR(logR + l) 



which, in view of the fact that xr = 1 in Ir-i and g is positive outside M x [—1,1] and 
bounded below by a constant otherwise, yields that for every R large enough, 



(4.6) 



/ 

JIr-1 



dg < dR(\ogR+ 1). 



This in turn implies — using (4.1 ) and the fact that ^\E\ 2 = g outside Rx [—1, 1] — the first 
(unsufficient) control 

/ |£| 2 < Ci_R(logi? + l). 

J{(x,y)\(x,0)<tUB P } 

Since the sum of the radii of the balls in B p intersected with any segment of M of length 1 
bounded by p < 1/8, we deduce by a mean value argument with respect to the variable x 
that there exists t £ [0, 1] such that 



(4.7) 



* (-f -t,v 



+ 



dy < CiR{logR+l). 
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Using now a mean value argument with respect to y, we deduce from (4.6) the existence of 
y R E [1,1 + V~R] such that 



(4.8) 



\E(x, y R )\ 2 + \E{x, -y R )\ z dx < C 1 VR(logR + 1). 



Next, we integrate div E = 2tt{v 



on the square 



R 

2 



using the symmetry property of Corollary 1.5 that 



"(1 





rVR 




R + 2t 


* J 


E 









R 
2~ 



t,y 



+ 



dy + 



f+t 



Vr,Vr]- We find 
E(x,y R )\ dx. 



Using the Cauchy-Schwarz inequality and (|4.7|)-(|4.8|), this leads for R large enough to 
(4.9) K^-i) 



R 



< 2 + C l R*l^\ogR + l + C l ^y R ^R{\ogR+l) < C l R 3 l\\ogR + 1), 
and then — since v{Ir] > K-^R-t/2) — to 

v(Ir) ~R> C 1 R 3/i log R. 

To prove the same upper bound for R — v(Ir) we proceed in the same way, but using a 
mean value argument to find some t G (—1, 0) instead of (0, 1) such that (4.7) holds, and then 
(4.8) also. We deduce as above that (4.9) holds and conclude by noting that since t S (—1,0) 
we have v{Ir) < v(I R _ t i^). This establishes (4.3). 

We may bootstrap this information: Indeed (4.3) implies in particular that v{I R ) — 
v{Ir-i) < CxR?/ A \ogR and thus we deduce from (O) that d45l) holds: 



W(E, XR ) 



XRdg 



< CiR 3/A log 2 R. 



Then since W(E,xr)/R — > W(E) as R — > oo uniformly w.r.t. E £ G and since g is both 
bounded from below by a universal constant and equal to ^l-E 1 ] 2 outside Ex [— 1, 1], we deduce 
(4.4), for R large enough depending on G. 



□ 



Definition 4.3. Assume v n = Y27=l^%i- Letting v' n = Y^i=x be the measure in blown-up 
coordinates and E Vn = —\/H' n , where H' n is defined by (1.21), we denote by g Vn the result of 



applying Proposition 4-1 to (u' n ,E Un ). 

4.2 Some preliminary construction lemmas 

The following lemmas serve to estimate the energy of explicit vector fields on boxes, which 
will be later combined to make test configurations. They are adaptations of |SSlj and rely 
on elliptic estimates. 

Lemma 4.4. Let K be the square f] 2 . Let (p G L 2 (dlC) and a £ L°°([— \, \\) be such 

L 

that J 2 L a(x) dx = — f dfC tp. Then, ao being the average of a over [— j, f], the solution (well 
defined up to an additive constant) to 



(4.10) 



-Ait = a<5f 
du 



Oi 



<4> 



in fC 
on dlC. 
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satisfies for every q £ [1,4] 



\S/u\ q < Cq ( a Q L 2 + L 2 \\a - a \ 



Proof. We write the solution u of (|4.10|) as u = u\ + U2 + «3 where 
(4.11) 



-Ami = ao^R in /C 
du i 



0i 



on <9/C 



where 99 is equal to on the vertical sides of the square and to #^ on both horizontal sides; 



(4.12) 



and 



-A112 = (a — oo)5r in K 
du2 



dv 



' -Au 3 = 
dv 



on dK 

in K 

if — (p on dK. 



The solution of (4.11) is (up to a constant) u\{x, y) = ^-\y\- Hence 



(4.13) 



I |Vm 

JK 



iq= (rmy L 2_ 



\ 2 J 



For U2, we observe that || (a— oo)^r[|w- 1 ««(x:) i s controlled, for any g < 00, by || a — ao || x°° f [ £. k]y 

Therefore, using elliptic regularity for (4.12), ([Vi^l I £«(*;) is controlled by \\a — ao\\ Loo nL 
and a scaling argument shows that for any q < 00 



(4-14) J \Vu 2 \ q < C q L 2 \\a - «o||^ ([ _. ^ 

Finally, in the proof of Lemma 4.16 of [SSlj it is shown that for any q £ [1,4] 



(4.15) 



IVu^ <C q L 2 -2M q L2{dlcy 



Combining (4.13), (4.14) and (4.15), we obtain the result 



□ 

Lemma 4.5. Let m be a positive constant and let K be a square of center 0, and sidelength 
1/m. Then the solution to 



-Af = 2ir(5 - m5$, 
dv 



in K 
on dK 



satisfies 



(4.16) 



lim 

r)->0 



>C\B(0,ri) 



|V/| 2 + 2vrlog77 



C — 7r log m 
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where C is universal, and for every 1 < q < 2 



(4.17) 



where C q depends only on 



[ \Vf\ q <C q mi-\ 
JK 



Proof. By scaling we can reduce to the case of m = 1. Then, it suffices to observe that 
f(z) = - log \z\ + S(z) with S e W X '°°(1C) and scale back. □ 



We note that W as defined in (1.6) still makes sense for vector fields satisfying div E = 
2tt(Y1 $x.i ~ m ) as long as curli? = in L)B(xi,r]o) for some t]q > 0. This is the notion we will 
use repeatedly below. 

Lemma 4.6. Let E\ and E<i be two vector fields defined in a rectangle 1Z of the plane which 
is symmetric with respect to the real axis, and satisfying 



(4.18) 
(4.19) 



div Ei = 2tt(^2 $xi - ait 

i 

div E2 = a2<5]R in 1Z 



in 1Z 



and curl E\ and curl E2 vanish near the Xi 's, for some distinct points x% G M and some bounded 
functions on the real line, a± and ai- Then, for q < 2 and q' its conjuguate exponent, we have 
Ei e L*(K) and E 2 G {K) and 



1 



W(Ei + E 2 , l n ) < W(Ei, l n ) + ||^i||L,(7t)||S2|| £ ,' (7e) + -WMh 



2 W 



where W is still defined as in (1.6). 
Proof We have 



/ 



\Ei + E 2 



lK\U i B(xi,i 1 ) 
By Holder's inequality we have 



L 



I -Ei 1 2 + I E2 1 2 + 2Ei ■ E2. 



L 



Ei ■ E2 



K\UiB(xi,ri) 



< ll-^i 1 1 iff CTe.) 1 1 -^2 1 1 (7^) 



The result easily follows. 



□ 



4.3 Proof of the screening result Proposition 3.4 



First we note that, in view of (3.10), if we assume G C Ai satisfies the hypothesis of Propo- 
sition 3.4 and < e < 1, then there exists yo > and Rq > such that for all E £ G, we 
have 



(4.20) VR> R , 

This motivates the following: 



\E\ 2 <e w R. 



inx{\y\>yo} 
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Lemma 4.7. Let E £ A\ satisfy (4.20), where < e < 1. Then, denoting z = (x,y), if 
\y\ > max(2yo, Rq), we have 



\E\ 2 (z) < C 



e 10 (\x\ + \y\ 



where C is universal. 



Proof. Each of the coordinates of E is harmonic in the half plane = {y > 0} since 
div E = cmlE = there. Therefore \E\ 2 is sub-harmonic. Thus, if B(z, \y\/2) C then by 
the maximum principle we have 



\E\ 2 (z 



Je 



\E\ 



B(*,|vl/2) 



Ify>2y , then B(z, \y\/2) C [x 



J|i,x + M] x [y 0) +oo)c 



1 2/1 , kl + 1 2/1 ] x [yo,+oo). 



Thus in view of (4.20), if |ac| + \y\ > Rq/2 we have 

\ E \\z) < i e ioN + H ! 

vr 

and the result follows, by symmetry with respect to the x-axis. 



□ 



Lemma 4.8. Let G satisfy the assumptions of Proposition 3.4 ■ Then for any E G G, any 
< e < I, 
such that 



< e < 1/2 and any i? /arge enough depending on G, e, we may find f G [| - ei?, =§ — |e-R] 



(4.21) 



where C depends only on G, and 



\E\ 2 < C£~ l 



(4.22) 



uniformly in G. 



({-t}u{t})x] 



1 



hm 7 -W(E,l I J = W(E) 

R^-oo At 



Proof. Take E £ G and apply Proposition 4.1 to E for some fixed < p < 1/8. We obtain 



a density g and balls Now, using (4.5) in Lemma 4.2 together with the bound (3.9), we 



deduce that if R is large enough depending on G then for any E £ G and denoting g the 
result of applying Proposition |4.1| to E we have 



(4.23) 



R/2-eR/2 p p 

/ (g(x,y)+g(-x,y)) dxdy<CR+ / dg<CR, 

x=R/2-eR JR JI r 



where C depends only on G and we have used the fact that g > —C everywhere and g > 
on the set {|y| > 1}. Then by using the fact that the radii of the balls in B p which intersect 
any given interval of length 1 is bounded by 1/8 we deduce that if R is also large enough 
depending on e, the measure of the set A of x G [R/2 — eR, R/2 — eR/2] such that {x, —x] x E 



does not intersect B p is bounded below by eR/4. This and (4.23) implies that the set T of 



t £ A such that f R (g(t,y) + g(—t,y)) dy < C/e has measure at least eR/8 if C is chosen 
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large enough depending on G, and (4.1 ) and the fact that g = ^\E\ 2 outside R x [—1, 1] imply- 
that (lOil holds for t G T. Thus 



(4.24) 



\{t e [12/2 - e!2, 12/2 - eR/2] | (4.21 ) holds}| > 



eR 



For (4.22) we argue as in [SSlj . Lemma 4.14. We let \ '■ [0, +oo) — > R be a monotonic 
function with compact support and let x(x,y) = xd^D- First we note that for any Radon 
measure p in M 2 we have 



/r+co 
X dp = ~ \ 
Jt=o 



X '(t)p(I 2 t) dt 



x'(u/2)p(I u )du. 



u=0 



This implies straightforwardly using the definition of W(E, x) that 

r+oo r 
2 / (W(E, ljj - 9 (/ u )) x'(u/2)di= / X^-W(^X)- 



On the other hand, by (4.2) and applying Lemma 4.2 (4.3), if \' is supported in [x,y] C 
[12/2,12], then the right-hand side is bounded by C(\x — y\ + 12 3 / 4 log 2 12) [|x'[|oc f° r an Y R 
large enough depending on G. Given now any p : M + — > R supported in [x, y] C [12/2, 12] we 
may consider the positive and negative parts p + and p_ , and and their primitives x+ an d X- 
with compact support, which are monotonic. Applying the above to x+ an d X- we fi n d 

/•+oo 

/ ljj - 5 (l u )) ^du^Cdx-yl+^^^HplU. 

iu=0 

Since this is true for any p supported in [12 — 2el2, 12 — eR], it follows by duality that 
(4.25) f \W(E, 1 7 J - g(I u )\ du < C(\x - y\ + 12 3 / 4 log 2 12). 

Then we divide [12 — 2el2, 12 — eR] into, say, [vR] intervals li, . . . , so that their length is 

equivalent to ey/R for large 12. Then, for large enough 12, on each such interval (4.25) implies 
that 

\W(E, ljj - g(I u )\ du < C12 3 / 4 log 2 12. 

Therefore the set of u G I k such that \W{E,lj ) - g(I u )\ < 8C12 3 / 4 log 2 12 has measure 
at least 7| lfc|/8 and since this is true on each interval, and changing variables, the set of 
t e [12/2 - el2,12/2 - eR/2] such t hat \W(E,1 T J - g{I t )\ < 8C12 3 / 4 log 2 12 has measure at 
least 7el2/16. Together wit h Q4.24| ), this implies the existence of t £ [12/2 - eR, R/2 - eR/2] 
such that both p~2~I| ) and ( |4.22[ ) hold. 

□ 



We now prove Proposition 3.4. Let G satisfy its hypothesis and choose < e < 1, and 
E € G. Applying Lemma 4.8 we find that if 12 is large enough depending on G, e, then there 
exists t £ [R/2 - eR,R/2 - eR/2] such that ( |4.21[ ) and ( |4.22[ ) hold. For any such integer 
12 G N we may also choose vr such that 



(4.26) 



e 3 R <y R < e 5/2 R. 
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Finally we choose s > t such that s — t £ [y R , yR + 1] and | - s 6 N, and start constructing 
the vector field Er. 



- Step 1: splitting the strip. We split the strip Ir 
strips (see figure below): let 

D = [-t,t] x [-y R ,y R ] 



R R] 
2 ' 2 J 



x M into several rectangles and 



D + = [t,s] x [-y R ,y R ] 
D- = [s,-t] x [-y R ,y R ] 
D+ = [a, R/2] x M 

= [-i?/2, -s] x R 
L»i = [-s, s] x ([y fl , y fl + R] U (-fl - -y R ]) 
Doo = [s, s] x ([R + j/b, +oo) U (-oo, -R - y R \) . 




First we let Er = E in Deb -^ii = in and below we are going to define Er on each 
of the other sets. 

Recall that from Corollary 1.5 we have that E(x,y) is the reflection of E(x,—y) with 
respect to the line {y = 0}. We denote by the trace E ■ V on the right-hand side of Dq 
where V is the outward-pointing normal to Dq, ip_ the same on the left-hand side, <ph the 
trace on the upper side of Dq (which by symmetry of the problem with respect to the real 



axis is equal to that on the lower side). From (4.21) and the Cauchy-Schwarz inequality, we 
have 



(4.27) 



\<P- 



+ \<p+r < 



c 



f-\ + / b+| < c 




and from Lemma 4.7 and (4.26) we have 

|2 ^ Ce 10 R(R + y R ) 



(4.28) 



■t,t] 



<Ce 4 



-t.t] 



\<Ph\< CeWR. 
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In addition, integrating the relation div E = 27r(^ peA 5 P — <5r) over Dq gives 
(4.29) i- H ^+ + J <p- + 2 J <pi\ = u([-t, t}) - 2t. 



- Step 2: defining Er in D + . We first define tpo := En ■ V on the boundary dD + , where V is 
the outward normal to D + . For a certain constant ipf to be chosen later we let 



(4.30) 



—(f+ ondD + ndDo, 
ifo = < on dD + n {x = s}, 

^+ on dD + n {y = ±y R }. 



Then, inside D + , we let E R = E\ + E2, where 



(4.31) 



div Ei = 2iriy~] 5 Xi - m + 5 R in P H 



Vl=l 



E!-v = a 



on dDj 



(4.32) 



div E2 = 2ir(m + — 1)5r in T> + 



E 2 -u = cp 



on dD + . 



Here, n + is an integer and m+ a real number which are defined by 



(4.33) 



n+ 



1 

2^ 



m + 



8-t' 



and for 1 < i < n + we have let 



8-t f. 1 

a?i = t H z + - 

n+ V 2 



Note that the above equations do not yield a uniquely defined E\ and E 2 - For (4.31 ) to make 



sense we need n + > while for (4.32) to have a solution we need to have 



(4.34) 



2vr (n+ - (s - t)) = / ^ = 2(a - 



which we take as the definition of 99^. The fact that n + > follows for f? large enough 
depending on e from the fact that s — t > e 3 R and (4.27), (4.28). 



- Step 3: Estimating the energy of Er in D + , D-. To compute the renormalized energy 
W(Er, 1d+) we need to define E\ and E2 more precisely. For E\ let us consider n+ identical 
squares {Ki}™^ with sidelength ^ = — , sides parallel to the axes and such that K{ is 
centered at X{. We define E\ restricted to Ki by applying Lemma 4.5 with m = m + and 
taking E% = —V/, while outside Uj-fQ we let E\ = 0. Since from Lemma 4.5 we have 
E\\ K , ■ v = on dKi, it holds that div E\ = ^ div E\\ K . and therefore (4.31 ) is satisfied by 
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Ei. On the other hand, still from Lemma 
(4.17) on the n + rectangles 



4.5 



we obtain by summing the bounds (4.16) and 



(4.35) 

and 
(4.36) 



lim 

j;->0 



D + \UiB(xi, V ) 



|£i| 2 + 27rlogr/ 



< n + (C — it log m + ) 



VI 



< q < 2, / |£7i|9 < C q n + . 



We define E2 by applying Lemma 4.4 in D + , hence with L = s — t, with the boundary 
data ifo and constant weight m = 27r(m + — 1). From (4.30) and (4.34) the hypothesis 
j m(x) dx = — j dD+ 93 is satisfied and applying the lemma yields 



(4.37) 



V2 



<q<4, [ \E 2 \ q <C q (\m + -l\ q {s-t) 2 + {s-t] 
Jd+ v 



2-2 

2 



V°IIl 2 (9D+) J • 



Using Lemma 4.6 we have, recalling that Er := E\ + E2 in L>+ and using (4.35), (4.36), 
(147371): 



(4.38) W(£ K , 1d+) < Cn + + C^ 1 ^ (\m+ - l|(s - i) 2 ^' + (a - t) 



\2/g'-l/2i 



+ C ( |m+ - l| 2 (s - t) 2 + (a - t)||^o" 2 



( Po\\L 2 (dD + ) 
h(dD+) 



for any 1 < q < 2 such that the conjugate exponent q' is less then 4. Now, from (4.33), we 
have using (4.27), (4.28), (4.26) and the fact that yn < s — t < yn + 1 we deduce that 

1 



(4.39) 



\n+-(s-t)\ < C e 2 ^ 




\m + - 1| < C 



1 



eVR 



+ 



and thus for R large enough depending on e, since s — t ~ y# for large i? and using (4.26) 
again as well as e < 1, 



(4.40) 



n + < Ce b ' 2 R, 



11 < 



C 



Moreover, from (|433|), (|4.34|) and (|4.28|) we find 
(4.41) |^+| < 



1 



2(s-t) 

Then, in view of ( |4.30| ), ( |4.27| ), 



2vr+ / ^ 



2//? 



eVR 



\ c Po\\l2(8d + ) 



< 



C 



Now we fix for instance g = 3/2, so that (/ = 3 and combining the above with (4.38), (4.40) 
we find that for R large enough depending on e, and denoting by C £ a positive constant 
depending on e but independent of R, 

W(E R , 1 D ) < Ce 5/2 i? + C e i? 2/3 x i? 2 / 3 " 1 ^ + <? £ 3/2 R 
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Thus for R large enough depending on e we find that 
(4.42) W(E R ,1 D+ ) <Ce z ' 2 R. 

An almost symmetric computation yields the same bound for W(Er, 1_d_). It suffices to let 

1 



(4.43) 



2(s - t) 



2?r 



ip + + / Lf- + 2 / Lp h 



and carry on the proof with minuses instead of pluses. The fact that n_ is an integer follows 



from the identity (4.29), the fact that 2s is an integer and the fact that u([—t,t]) € N. 

(s-t) 



Moreover the definition of n_ implies that 
{a-t)- — ( [ <p-+ [ ip h 



71- 



or n_ 



1 

2^ 



<P-+ <Ph 



+ 1, 



hence n_ is positive if R is large enough and (4.39) holds for ra_ as well. The rest of the proof 
is unchanged. 

- Step 4 : defining Er over D\. We need only consider the intersection of D\ with the 
upper half-plane (and then extend by reflection). We let tpo be equal to — tp^, —iph, ~ ( Pf l i 
respectively, on the intersection of dD\ with dD_ , 8Dq , dD + , respectively. On the remaining 



three sides of dD\ we let fo = 0. From (4.34) and its equivalent for ra_ and the fact that 
n± = (s — t)m± we have 



- J fo = vr(n+ + n_ - 2(s - t)) + ^ J + \ j f-+ J 



<Ph, 



and then (4.33), (4.43) imply that the integral of cpo is zero. 



Thus there exists a harmonic function u in D\ with normal derivative (po on dD\, we 



let Er = Vu on D\, Using (4.41) — which holds for ip h as well — and (4.28) we have 
Ce 4 + Ce 4 R/y R hence 

(4-44) lko|||2 (9Dl) < C {e 4 + e 4 R/y R ) < Ce. 

Then standard elliptic estimates yield as in Lemma |4.4| that 



(4.45) 



/ |^r| 2 < CJJl^ollia < CRe, 



where we have concluded by (4.26). 



- Step 5: defining Er over D+. The construction will be entirely parallel in D e . We note 
that D+ is an infinite strip of width R/2 — s and we have chosen s so that this quantity is an 
integer. We can thus split this strip into exactly R/2 — s strips of width 1. On each of these 
strips we define Er to be equal to for \y\ > | and for \y\ < \ (i.e. in a square of sidelength 
1) we choose it to be V/ where / is given by Lemma 4.5 applied with m = 1. Since Er -v = 



on the boundary of each of these squares, no divergence is created at the interfaces, and the 
resulting Er satisfies div Er = 27r(^ pgA <5 p — <5k). In addition in view of (4.17) the cost in 
energy is equal to a constant times the number of strips, that is 



(4.46) 



W(E R , l De ) < C\R/2 -s\< CeR. 
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- Conclusion. We have now defined Er over the whole strip Ir. It satisfies items ii) and iii). 
The main point is again that as long as E ■ v is continuous across an interface it creates no 



singular divergence there. Combining (4.46) with (4.22), (4.42) and (4.45), En also satisfies 



(3.11). This concludes the proof of Proposition 3.4 



5 Proof of Proposition 3.1 



The construction consists of the following. First we select a finite set of vector fields Ji, . . . , Jn 
(N will depend on e) which will represent the probability P(x, E) with respect to its in- 
dependence, and whose renormalized energies are well-controlled. Since P is 7a(x) -invariant, 
we need it to be well-approximated by measures supported on the orbits of the Jj's under 
translations. It is also during this approximation process that we manage to select the Jj's as 
belonging to a part of the support of P of almost full measure for which the extra assumption 



(3.10) holds and the screening can be made. Secondly, we work in blown-up coordinates and 



split the region £' (of order n size) into many intervals, and then select the proportion of 
the intervals that corresponds to the weight that the orbit of each Jj carries in the approxi- 
mation of P. In these rectangles we paste a (translated) copy of J{ at the appropriate scale 
(approximating the density mo' by a piecewise constant one and controlling errors). 



To conclude the proof of Proposition [3Tj we collect all of the estimates on the constructed 
vector field to show that its energy w n is bounded above in terms of f W dP and that the 
probability measures associated to the construction have remained close to P. 



5.1 Abstract preliminaries 



We repeat here the definitions of distances that we used in [SS2]. First we choose distances 
which metrize the topologies of L^ C (IR 2 , M?) and B(X), the set of finite Borel measures on 
X = E x Lf oc (M 2 ,]R 2 ). For E U E 2 G Lf oc (M 2 , M 2 ) we let 



d p (E 1 ,E 2 ) = J2r 



k=l 



. \\Ei — E 2 \\Lp(B(o,k)) 
1 + 1 1 -El - -E , 2||lp(B(0,A:)) ' 



and on X we use the product of the Euclidean distance on £ and d p , which we denote dx- 
On B(X) we define a distance by choosing a sequence of bounded continuous functions {<fk}k 
which is dense in Cb(X) and we let, for any /ii,/U2 £ B(X), 



k=l 



-k K^fc^l ~ ^2)| 
1 + \ (<Pk,fJ>l ~ M2) 



where we have used the notation {ip,fj) = j ipdfi. 

We will use the following general facts, whose proofs are in |SS2} Sec. 7.1]. 

Lemma 5.1. For any e > there exists rjo > such that if P,Q G B(X) and \\P — Q\\ < rjo, 
then d(P, Q) < e. Here \\P — Q\\ denotes the total variation of the signed measure P — Q, i.e. 
the supremum of P — Q) over measurable functions ip such that \<p\ < 1. 



In particular, if P 
d B (P,Q)<e. 



Y^i a i S xi and Q = witn ~ &l < Vo, then 
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Lemma 5.2. Let K C X be compact. For any e > there exists rji > suc/i i/ia£ if 
x E K,y £ X and dx(x, y) < rj\ then ds(5 x ,5 y ) < e. 

Lemma 5.3. Let < e < 1. If n is a probability measure on a set A and f,g:A^X are 
measurable and such that ds(S^f x \,6 g f x \) < e for every x G A, then 

d B (f#^g#fi)<Ce(\loge\ + l) 

where # denotes the push-forward of a measure. 

The next lemma shows how to, given a translation-invariant probability measure P on 
L^ oc (IR 2 , IR 2 ), select a good subset G £ and vector fields J, of L^ oc (IR 2 , M 2 ) to approximate it. 
It is essentially borrowed from |SS2| except it contains in addition the argument that ensures 



that we may choose G £ to satisfy the assumption (3.10) needed for the screening. 



Lemma 5.4. Let P be a translation invariant measure on X such that, P-a.e., E is in A\ 
and i 
that 



and satifies W(E) < +oo. Then, for any e > there exists a compact G £ C L^ c (M 2 ,IR 2 ) such 



i) Letting < t/q be as in Lemma 5.1 we have 



(5.1) P(SxG e c )<min(r ?0 2 ,? ? oe). 
ii) The convergence ( |1.8[ ) is uniform with respect to E G G £ . 

Hi) Writing div E = 2tt(i/e — 1); both W(E) and ve{Ir)/R are bounded uniformly with 
respect to E £ G e and R > 1. 

iv) Uniformly with respect to E G G e we have 

(5.2) lim lim I / |£| 2 = 



Moreover, (5.1) implies that for any R £ > 1 there exists a compact subset H £ C G £ such that 
v) For every E G H £ , there exists T(E) C ImR E such that 

(5.3) \T(E)\ < CR £ rj and A <£ T(E) =>- 9\E G G £ . 

vi) We have 

(5.4) d B (P, P') < Ce(|log e\ + 1), where 

P '= I ~ T^Tr i / 6 X ® 5 dA dP(z, J5) 

m (x); flE \r(£) 



^ = / k®* a moix)E dP(x,E) 
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Finally, there exists a partition of H £ into U^iJ* satisfying diam (HI) < 773, 
such that 



where r/3 is 



(5.5) E G H £ , d p (E,E') < 773, m G (0,m], A G mI Re \T(E) 
and there exists for all i, Ei G suc/i i/iai 



dB{&<j m 6 x E,o~<j m 6 x E') < e; 



(5.6) 



< infVF + e. 



Proof. The lemma is almost identical to Lemma 7.6 in [SS2J, except for item iv). The proof 
in [SS2] is as follows: First one proves that there exists G £ satisfying items ii) and hi) with 
P(E x G £ c ) arbitrarily small, in particular one can choose it so that (5.1) is satisfied. Then 



one deduces from (5.1) the existence, for any R £ > 1, of a compact subset H £ C G £ satisfying 



the remaining properties. The only difference here is that we must check that there exists G £ 
with P(S x G £ c ) arbitrarily small satisfying not only items ii) and hi), but iv) as well. Then, 
the proof of the existence H £ C G £ satisfying the remaining properties is exactly as in [SS2]. 

Of course, by intersecting sets, it is equivalent to prove that ii), hi), and iv) can be satisfied 
simultaneously or separately, on a set of measure arbitrarily close to full. The proof in |SS2j 
shows that this is possible for ii) or hi), it remains to check it for iv). For this we consider 
G n = {E I W(E) < n}. Then G n is a translation-invariant set since W is a translation- 
invariant function, and therefore by the multiparameter ergodic theorem (as in |BeJ), and 
since P is translation- invariant, we have 
(5.7) 



where xr 



\E\ 2 dP(x,E) = 
-l,l]x{\y\>y } J JXxG n 

Then, using Lemma [472 



* lr 



defined in Proposition 



L1L 



4.1 




dP(x,E), 



xr\e\ 2 
-x{\y\>yo} j 

and using the fact that the g there was 



hence is equal to \\E\ 2 on E x {\y\ > 1} we deduce from (4.5) and 
the fact that g is bounded below by a constant independent of E that 

1 



lim 

R-^oo R 



Xr\E\ 2 < C(1 + W{E)) < Cn 

^{\y\>yo} 

holds for every E G G n with n > 1. 

It follows that for every fixed n > 1 the family of functions 



\E\ 



-l,l]x{\y\>y } 



decreases to on S x G n as yo — > +00, and is dominated by the bounded, hence P-integrable, 
function tp\. Lebesgue's theorem then implies that their integrals on E x G n converge to 0, 



hence in view of (5.7) that 



S x G n 




x{\y\>yo} 



XR \E\ 2 dP(x,E), 



+00. Fatou's lemma then implies that (5.2) holds for P-almost every 



tends to as yo 
(x,E) G £ x G n . 

Since W(E) < +00 holds for P-a.e. (x,E), we know that P(S x G n ) — > 1 as n — > +00 
therefore the measure of S x G n can be made arbitrarily close to 1, and then Egoroff's 



theorem implies that by restricting G n we can in addition require the convergence in (5.2) to 
be uniform. □ 



30 



5.2 Construction 



In what follows £' = nS, mo'(x) = mo(x/n): we work in blown- up coordinates. In view of 
assumption ( |1.16| ), we may assume without loss of generality that £ is made of one closed 
interval [a, b] (it is then immediate to generalize the construction to the case of a finite union 
of intervals). In that case £' = [na,nb]. Let m > be a small parameter. For any integer 
n we choose real numbers a n and b n (depending on m) as follows: Let a n be the smallest 
number and b n the largest such that 

nm nm 

(5.8) a n >na-\ b n < nb 

(5.9) / m '(i)rfieN 

/•nfe 

(5.10) / m '(x)dx£N 



bn 



(5.11) / m '(x)&G# 



fin 



where g e is an integer, to be chosen later, and 7 is the constant in (1.17). By (5.8) and 



assumption (1.16), we are sure to have mo' > m in T,' m := [a n , b n ]. This fact also ensures that 



nm 1 nm q £ 

(5.12) \a n -na\<^ r -\ \b n — nb\ < — H . 

7 Z m 7 UL 



We also denote := ^' m . 

Let P be a probability on £ x L^ c (IR 2 , M 2 ) which is as in the statement of Proposition 
3.1 Our goal is to construct a vector field E n whose W energy is close to j W dP and such 
that the associated P n (defined as the push-forward of the normalized Lebesgue measure on 
£ by x 1— > (x, E(nx + •))) well approximates P, a probability measure as in the statement of 
Proposition |3.1| 

In [na,a n ] and [b n ,nb], we approximate mo{x)dx by a sum of Dirac masses at points 
appropriately spaced, and build an associated E n , whose contribution to the energy will 
shown to be negligible as m — > 0. We leave this part for the end. 

For now we turn to [a n , b n ], where we will do a more sophisticated construction, approach- 



ing P via Lemma 5.4 and using Proposition 3.4 The idea of the construction is to split the 



interval [a n , b n ] into intervals of width ~ q £ Re, where q e is an integer and R £ a number, both 
chosen large enough, and then paste in each of these intervals a large number of copies of the 



(rescaled) truncations of the Jj's provided by Proposition 3.4, in a proportion following that 
of P. 

-Step 1: Reduction to a density bounded below. We have 

P = j 5 X ® 5 amo{x)E dQ(x,E), where Q = J S x (g> 6 Ul/mo(x)E dP(x, E). 

Moreover, since the first marginal of P is the normalized Lebesgue measure on S and since 
|S m | ~ |S| as m — > 0, we have 

(5.13) lim d B (P,P) =0, 

m— >0 
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where P is defined by 



(5.14) P= / 6 x ®6 amo(x)E dP(x,E), with P 



S m xLf oc ( 



5 x ®5 ai/mo{x)E dP(x,E). 



Clearly P is is Txr x ) -invariant since P is, and in particular it is translation-invariant. In 
addition, for P-a.e. (x,E), we have mo(x) G [m>^]) a situation similar to [SS2] where the 
density was assumed to be bounded below. 



Step 2: Choice of the parameters. Let < e < 1. We apply Lemma 5.4[ w hich provides us 

2 ,IR 2 ). We may then apply Proposition 



with a compact set G £ C L p 



loc\ 



3.4 



with the same e 



to G £ . It provides us with Rq > such that for any integer R > Rq, and roughly speaking, 
any E £ G £ may be truncated on Ir with an e-cost in energy. Applying Lemma 5.2 on the 
compact set {o- m E : m G [m, m], E 1 G Cr £ }, there exists r/i > such that 
(5.15) 

m G [m,m], E G G £ ,E' G Lf oc (R 2 ,IR 2 ) and d p (E, E') < r/i =► d B (^ m £,^ m B') < £• 
Then we define i? e to be such that mi? £ > i?o and such that for any E, E' G L^ C (M 2 , M 2 ), 
(5.16) E = E' on 7„ ei?e => d^, £') < Vl . 



Going back to Lemma 5.4, we deduce the existence of H £ C G £ , of N £ G N and of 
{Ei} x <i< Ne satisfying Q, Q, @ and @. 
Finally, we choose g e G N sufficiently large so that 



(5.17) 



N, 



N. 



» x max W(a m Ei) < e. 

2 0<i<N E 
m£ [m,m] 



- Step 3: construction in [a n , b n ]. We start by splitting this interval into subintervals with 
integer "charge". This is done by induction by letting to = a n and, t^ being given, letting 
tfc+i be the smallest t > tk + q £ R £ such that i mo'(x)dx G g £ N. By (5.11) there exists 
K G N such that tx = b n , and 



(5.18) 



K < 



< 



n(b — a) 



q £ R £ q £ R £ 

Since mo' > m in [a n ,6 n ], it is clear that ty. — (q £ R £ + ifc-i) < q £ wT l "• To summarize and 
letting Jfc = [tfc_i,tfc], w e thus have 



(5.19) 



\h\ £ [qeRe,qe{Re +m / m '(z) G q £ N. 



In each we will paste & copies of a rescaled version of E^ , where 



q £ (b n - a n ) 



\h\ 



~Pi,k 



Pi,k = P ( -h x H\ 
n 



[•] denoting the integer part of a number. Because the first marginal of P is the normalized 
Lebesgue measure on and since [a n , b n ] C C [a, b], and Ufc/fc = K, bn], we have that 



141 



\h\ 



n[b -a) f-f 6 n - a„ 



i=i 
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and therefore Yli=i n i,k — Qe- Also, using in particular (5.17), 

141 



'^i,k Pi,k 



N £ 

< — < Vo- 
le 



We divide Ik into q £ subintervals with disjoint interiors, all having the same width E 
[R £ , R £ +WT 1 ]. Then for each 1 < i < N e we let Ii t k denote a family consisting of nj & of 



these intervals. This doesn't necessarily exhaust Ik since n i,k — Qs so we let n o,k = 

We define to be the average of mo' over From (5.19) we have mk\Ik\ E hence 
for each I E I^k we have := |m&I| E N, and R E [m,kRe,mk(R £ + HI )]■ We then apply 
Proposition 3.4 in i# to the vector field Ei, which yields a "truncated" vector field En defined 
in Ir, where R = \mi~I\. If I E Xo,& we apply the same procedure with an arbitrary current 
Eq E A\ fixed with respect to all the parameters of the construction. 

We then set 

E n \ x ) = Vl/m h Ei,l{xi + •) 

on each interval I E Zj where xj is the center of /. The next step is to rectify the weight in 
En \ For this we let IZk be the square Ik X (— |ifc|/2, |Ifc|/2) and let Hk be the solution to 

-AHk = 27r(mo' — m^) inT^fc 

dH k 



dv 







on dTZk- 



From Lemma 4.4 applied with tp and mo equal to zero, and using the fact that mo is assumed 
to belong to C°2 , we have for any q E [1,4], 



(5.20) 
We then define 



/ \VH k \ q < C g |/ fc | 2 ||mo' - m k \\ q Loo{Ik) < C g |4| 2 ||m | 



£70,1 



_ 2 

n 2 . 



M 2 ) 



VF fe in Tlfe 
^ in JfcVT^fe 

E n = E^+E^ in 4. 



Using Lemma 4.6 and (5.20) we deduce using (5.19) that 

(5.21) W(E n , l fk ) < W{E%\ l f J + o n (l), as n -> oo, 

where o n (l) tends to zero as n — )• oo and depends on e,m > but not the interval Ik we 
are considering. Summing (5.20) for 1 < k < K and in view of (5.18) we find that for any 
9 €[1,4] 



(5.22) 



f {Eg* - E n \« < Cq^rnfl 1 - 

J[a n ,b„)xR 



(5.23) 



On the other hand, in view of the construction and the result of Proposition 3.4 we have 

/ N £ \ 



W(Ep,l Ik )<\I k \lJ2 



,i=0 



-W(a mK Ji) + Ce 
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Then, following the exact same arguments as in [SS2, Sec. 7] which we do not reproduce 
here (the only difference is that the rescaling factors \fn there should be replaced by n), 
thanks to ( |5.16 )— ( 5.15 )— ( 5.17 ) we find that we can choose C\ in (5.17) such that 



(5.24) 
where 



d B (P, P')<C 'e(|loge| + l) 



A 



P' 



w k= i 



h 



and sta nds for P^ in |SS2|. Sec. 7]. Also, and again as in |SS2j . since (5.22) holds, and from 
Lemma 



5.2 



we may replace E n with E n at a negligible cost, more precisely for any large 
enough n we have 



(5.25) 
where 



d B (P,P")<Ce(\loge\ + l) 



K 



P" 



k =l 



- Step 3: construction in [b n ,nb]. The construction in [na, a n ] is exactly the same hence will 
be omitted. We claim that there exists E n defined in [b n ,nb] x IR such that 



div En 



(5.26) 



2vr(^^ 



m 



E n ■ V 



in [b n , nb] x R 

on d([b n , nb] x R) 

and 

(5.27) W(E n , l [6ni „ 6]xR ) <Cn(m + o n (l)) , 

where C may depend on j,fn and e. To prove this claim, let so = b n and for every I > 1, let 
si be the smallest s > s/_i such that mo'(x) dx = 1. Since (5.10) holds, this terminates 

at some sl = nb with L = < m\nb — b n \. We then set xi to be the middle of [s/_i,s;]. 

We let ui be the solution in the square IZi := [s/_i, si] x [—^(si — s;_i), ^(si — si^i)] 



-Aui = 2tt(8 Xi 
Qui 
dv 



m 







in IZi 
on dlZt- 



This equation is solvable since, by construction of the s/'s, the right-hand side has zero 
integral. Then for each / we let E n = —Vui in 7Z[, and let E n = in [b n ,nb] x R \ U[7Zi. 
Clearly E n satisfies (5.26). 

To estimate the energy of u\ we let u\ = Vi + wi where, letting mi = 4j S; i s ^ too', 

miSyi) in IZi 
on dTZi, 



-Avi = 2tt(5 x 
dvi 
dv 
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-Awi = 2Tr(mi - m ')4 in Hi 
dwi 



dv 







on &R-i, 



From Lemma 4.4 and Lemma 4.5 we find, choosing for instance q = 4 so that 5 £ [1)4] and 
q'<2, 



n, 



iV^I^C^-^ifllmz-mo' 11 '' 



L°°([aj_i,s,]) 



and 



W(«,,l Wi ) = C7-?rlogmj, / |V^| 9 ' < Cm, 9 '- 2 . 



From (1.17) and Lemma 4.6, since £7 n = — (Vu; + Vw;) in 7?./, we have 



(5.28) 1^(^,1^) <C- tt log mj + C 

+ C||m i -m , ||^ oo([s; _ iiS;]) (ai-ai-ij 



Using (1.18), 



I 'II s n\\ 11 s«-s;-i 2 

m; — mn roo / r n < C mn n l 



Replacing in (5.28) and letting q = 4 we deduce that 



W%E„, W < C-7T log m, + C 



- sz-i (gj - g;-i) 3 



Then, summing with respect to / — using the fact that from (5.12) we have ^ ~ s i\ < 
Cnm(l + o n (l)), the fact that the integral over s;] of mo' is 1 and that from (1.17) we 
have (sj — sj_i) < n5 — we find 

W(E n , l[b n>n b]xM.) <C^J^ m '(x) - m '(x) logm '{x) dx + no„(l)^ < Cn(m + o n (l)), 



since mo' — mo' log mo' is bounded by a constant depending only on mo and using (5.12). 
This proves (5.27) 



Step 4-' conclusion. Once the construction of E n is completed, the proof of Proposition 3.1 



essentially identical to that of [SS2, Proposition 4.1], which is its 2-dimensional equivalent, 
except that the scaling factor yfn there must be replaced by n. We only sketch the proof 
below and refer to the specific part of |SS2| for the details. 

The test vector-field E n has now been defined on all [na, nb] x M. It is extended by 
outside, and is easily seen to satisfy the relation div E n = 2ir{u' n — mo') for u' n = Yli=i $x'i a 
sum of Dirac masses on the real line. Combining (5.27) with (5.21), ( 5.23[ ) and (5.6), we have 



N e 



W(E n , 1 R2 ) < 141 ( Yl —W(a mK Ei) + Ce + o n (l) ) + Cn(m + o n (l)). 



,i=0 
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Letting n — > oo and then m — > 0, we see that the error term on the right-hand side can 
be made arbitrarily small, say smaller than Ce. On the other hand, the reasoning of [SS2J, 
Step 2 in Paragraph 7.4, shows that 



Yl I 1 *! S —W{<Jm K Ei) < |E'| / W(£7) dP(x, E) + Cn(e + o n (l)) , 
fc Vi=o 9e ) 3 

so that taking n larger if necessary we obtain 

(5.29) ^-W(E n ,l R2 ) < j W(E)dP(x,E) + Ce. 

Then arguing as in Paragraph 7.4, Step 3 of |SS2j we obtain that, letting (xi, • • • , x n ) be 
the rescalings to the original scale of the points x\ i.e. Xi = x[/n, we have for n large enough 

1 |S| f 

limsup — (w n (xi, . . . ,x n ) — n 2 J r (no) + nlogn) < — / W(E) dP(x,E) + Ce. 

n— >oo n IT J 

Also letting P n be the push-forward of r^idx |s by the map x \-t (x, E n (nx + •)), it is easy 



to see that d&(P" , P n ) < Cm. In view of (5.13) and (5.25), and taking m small enough, for 
any given e > 0, we can achieve 

d B (P,P n )<Ce. 



This proves that items i) and ii) of Proposition 3.1 are satisfied by (xi, . . . , x n ) and E n . 



Then, the perturbation argument of Paragraph 7.4, Step 4 in |SS2| shows that there exists 
5 > and for each n a subset A n C W 1 such that \A n \ > n\(5/n) n and such that for every 



(Vl, ■ ■ ■ iVn) G A n there exists a corresponding E n satisfying (|3.l| and (|3.2[) . This concludes 



the proof of Proposition 3.1 



6 Proofs of the remaining theorems 
6.1 Preliminary bounds on Z@ 



The crucial fact is that there are good lower bounds for F n . This follows from Proposition 4.1 



Even though we will not use the following result in the sequel, we state it to show how we 
can quickly derive a first upper bound on Zn from what precedes. 



Proposition 6.1. We have 
(6.1) 
and 
(6.2) 

where o(l) 



logK? < Cn(3 + n(log |S| +o(l)) 

logZf < -Pn 2 T(n ) +/3nlogn + C«£ + n(log|£| +o(l)) 
as n —> oo uniformly with respect to (3 > /3q > 0, and C depends only on V . 



The proof uses two lemmas. 
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Lemma 6.2. For any v n = YH=i $xi> we have 



(6.3) 



F n iy n ) 



l 



dgv n 



where F n is as in (1.23) and g Vn is the result of applying Proposition 4-1 



to v n . 



Proof. This follows from (4.2) applied to XRi where xr is as i n (1-7). If R is large enough 



then #{p € Supp(^) | B(p, C) D Supp(Vx) / 0} = and therefore (4.2) reads 



W(E Un , XR ) 



XRdg„ n . 



Letting R — > +oo yields W(E Un , 1 K 2) = J dg Un and the result, in view of (1.22). 
The following lemma has the same proof as in |SS2|. Lemma 3.5]. 



□ 



Lemma 6.3. Letting v n stand for Y^i=i we have, for any constant a > and uniformly 
w.r.t. > fa > 0, 



(6.4) 



lim 

n— >oo 



- a ? n ^ dUn dxi...dx r 



ISI. 



Proposition 6.1 is then proved exactl y as in |SS2j . Proposition 3.3. Note that a lower 
bound for Z% was obtained in Corollary 3.3 To complete the proof of Theorem [3] it then 
suffices to prove the corresponding upper bound. 



6.2 Lower bounds for F n and consequences 

We start with a result that shows how F n controls the fluctuation u n — n/xo- 

Lemma 6.4. Let v n = ^22=1^' ^or any interval L of width R (possibly depending on n) 
and any 1 < q < 2, we have 

1 

12 



\\v n - n/i ||vp-i.9(j) < C q (l + R)i 2 na \ F n (v n ) + 1 

Here W~ l ' q is the dual of the Sobolev space W^ q with 1/q+l/q' = 1. 
Proof. In |SS2t Lemma 5.1], we have the following statement 



11 1 



-l,9( Bfl ) <C q (l+R A )« 5 n5 (F n (v n ) + 1 



The proof is based on [SeTiJ which works in our one-dimensional context as well, thus the 
proof can be reproduced without change. It is immediate to deduce the result. □ 

We now turn to bounding from below F n . 

The proof is the same as in [SS2, Sec. 6], itself following the method of |SSlj based on 
the ergodic theorem. We just state the main ingredients. 

Let {v n }n and P Vn be as in the statement of Theorem |4j We need to prove that any 
subsequence of {P Vn }n has a convergent subsequence and that the limit P is a r A ( x )-invariant 
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probability measure such that P-almost every (x,E) is such that E G -4. mo ( x ) and (1.29) 
holds. Note that the fact that the first marginal of P is follows from the fact that, 

by definition, this is true of P Vn . 

We thus take a subsequence of {P Un } (which we don't relabel). We may assume that it 
has a subsequence, still denoted u n , which satisfies F n (y n ) < C, otherwise there is nothing to 
prove. This implies that v n is of the form Yl?=i 'Wn- We let E n denote the electric field and 

g n the measures associated to v n as in Definition |4.3l As usual, — y -j &nxi n • 

we have 



6.4 



A useful consequence of F n (u n ) < C is that, using the result of Lemma 

(6.5) —v n —?• /io on M.. 

n 

We then set up the framework of Section 6.1 in [SS2] . We let G = £ and X = A4+ x 
L^ oc (]R 2 , M 2 ) x Ai, where p £ (1, 2), where 7W+ denotes the set of positive Radon measures on 
R 2 and M. the set of those which are bounded below by the constant — Cy := — C(||mo||oo + 1) 
of Proposition 4.1, both equipped with the topology of weak convergence. 

For A G M and abusing notation we let 6\ denote both the translation x i— > x + A and the 
action 

O x (v,E,g) = (e x #v,Eoe x ,e x #g). 
Accordingly the action T n on S x X is defined for A G K by 

TZ(x,v,E,g) = (x + ~,ex#v,Eoe x ,9x#g 

Then we let x be a smooth nonnegative cut-off function with integral 1 and support in [—1, 1] 
and define 



-/ X(t)dg{t,s) if (v,E,g) = 6 nx (v n ',E n ,g n ), 
7T Jm.2 



(6.6) t n (x,u,E,g) = {n 

^ +oo otherwise. 

Finally we let, 

(6.7) F n (u,E,g) =j i n (x,9 xn (is,E,g)) dx. 

We have the following relation between F ra and F n , as n — > +oo (sec [SS2, Sec. 6]): 



< rkFn^n) + o(l) if (p, E, g) = (V, E n , g n 



|S| J 

+oo otherwise. 



(6.8) F n (v,E,g) is 



The hypotheses in Section 6.1 of [SS2] are satisfied and applying the abstract result, Theorem 
6 of [SS2], we conclude that letting Q n denote the push-forward of the normalized Lebesgue 
measure on X by the map x \- > (x, 9 nx (v n ', E n ,g n )), and Q = lim n Q n , we have 

(6.9) lim inf ^-Fn(v n )>- [ W(E) dQ(x, u, E, g) 

n |2j| TT J 

and, Q-a.e. (E,u) G An (a:)- 
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Now we let P n (resp. P) be the marginal of Q n (resp. Q) with respect to the variables 
(x, E). Then the first marginal of P is the normalized Lebesgue measure on E and P-a.e. we 
have E G A mo ( x y in particular 



W(E) > min W = mo(x) 



A 



m (x) 



min W 

Al 



7rlog mo(x) 



Integrating with respect to P and noting that since only x appears on the right-hand side 



we may replace P by its first marginal there, we find, in view of (1.27) that the lower bound 



(1.29) holds. 



6.3 Proof of Theorems H], [5] and [3J, completed 

As mentioned above, Part B of Theorem [4] is a direct consequence of Proposition 3.1 



see 



Corollary 3.2 



Part C follows from the comparison of Parts A and B: for minimizers, the chains of 
inequalities (1.29) and (1.30) are in fact equalities and that ^ j W dP must be minimized 



hence equal to a, cf. (1.35). Since (1.29) follows from (6.9), the inequality (6.9) must also be 
an equality and moreover we must have limn-s.00 J £ du n = 0. 

This completes the proof of Theorem [4j The proof of Theorems [5] and Theorem [3] is 
identical to |SS2j except again the replacement of the scaling y/n by n and \A n \ > n\(-K5 2 /n) n 
by \A n \ > n\(n5/n) n . 



6.4 Proof of Proposition 2.1 



The result of Proposition 2.1 is a consequence of Proposition 3.4 



First, applying Part C of Theorem |4j there exists a translation-invariant measure P on 
£ x L^ c (IR 2 , M 2 ) such that P-a.e. (x, E) is such that E minimizes W over A. mo t x ) ■ Then, taking 
the push-forward of P under (x,E) \— > ai/ mo ^E, we obtain a probability Q on L^ oc (M. 2 , M 2 ) 
such that Q-a.e. E minimizes W over A.%. 



Applying Lemma 5.4 to Q, we find that Q-a.e. E is such that E G A\, such that (3.10) 



to 



holds, and such that W(E) = min^ W. Choose on such Eq and apply Proposition 3.4 
G = {Eq}. We find that for any given e > and any integer R large enough depending on e, 
there exists Er defined on Ir such that Er ■ V = on din and W(Er, 1j r ) < R(W(E ) + e). 
This Er can be extended periodically by letting Er(x + kR, y) = Er(x, y) for any k G Z. 
From the condition Er ■ V = on 8Ir we find that, letting A C Ir be the locations of 



the Dirac masses in div Er, we have div Er = 
Moreover 

W(E r ,1 Ir 



ME 



W(E_ 



\Ir\ 



< W{Eq) + Ce. 



where A R = A + PZ. 



There remains to make Er curl-free. Following the proof of Corollary 4.4 of |SSlj we let 
Er = E r + V^-J'r in Ir where f R solves -Af R = cuvIEr in I R and f R = on dI R . Then, 
div Er = div P^ and curlPfj = in We can thus find Hr such that Er = VPr in Ir. 
It also satisfies VHr ■ v = Er ■ v = Er ■ V + V^Jr ■ v = on d Ir. We may then extend Hr to 
a periodic function by even reflection, and take the final Er to be VHr. This procedure can 
only decrease the energy (arguing again as in |SS11ISS2"] ): we have W(Er, 1j ) < W(Er, 1j ) 
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since 



/ \VH R -V L f R \ 2 - I \VH R \ 2 

Jl R \UB(p, V ) Jl R \UB(p, v ) 

= -2/ VH R -V ± f R + f \Vf R \ 2 . 

Ji R \uB(p, v ) Ji r \ub(p, v ) 

It can be checked that the last two terms on the right-hand side converge as 77 — > to the 
integrals over I R . Also integrating by parts, using the Jacobian structure and the boundary 
data, we have fj VH R • V f R = 0. Therefore, letting r\ — > in the above yields 

W(E r ,1i r )-W(VH r ,1 Ir ) >o. 

We deduce that W(E R ) < W(E R ) < min^ W + Ce, with E R a 2i?-periodic (with respect to 
the variable x) test vector field belonging to A±. The result follows by a standard diagonal 
argument. 



6.5 Consequences for deviations: proof of Theorem [6] 

The proof relies on the following proposition, whose proof is much shorter than in |SS2j , due 
to the simpler nature of the one-dimensional geometry. 



Proposition 6.5. Let u n = Y17=i^ x i' an( ^ 9 v n ^ e as * n Definition 4-3. For any R > 1, for 

any xq G M, denoting 



D(x ,R) = v n \ Br(x )) - nno [Br(x ) 



we have 



[ dg Vn > -CR + cD(x , R) 2 min (l, \^^A) 

Jb 2R (x' ) \ R J 

where c > and C depend only onV.^\ 

Proof. Two cases can happen: either D{xq, R) > or D(xq, R) < 0. 

We start with the first case. Let us choose r = min ( 2, 1 + ^( a o.^) \ j denote T 



{r G [R, tR], B r (x' ) fl Bp = 0}, where B p is as in Proposition 4.1 By construction of B p 
and since p < ^, we have \T\ ~> \(t — 1)R. We then follow the method of "integrating over 
circles" introduced in |SS3j : let C denote {x G B tR (x' )\B r (x' ), \x — x' Q \ ^ T}. 
For any r G T, since dB r (x' ) does not intersect Supp{v' n ), we have 

(6.10) / E Vn -v= I div E Vn = 2nu' n (B r (x , )) - [ m 

JdB r (x' ) JB r (x' ) JB r (x' ) 

> D(x , R) - 2(r - l)R\\m \\L°° > ^D(x , R) 
by assumption and by choice of r. Moreover, for any r G T, we have, by Cauchy-Schwarz, 



8B r (x' ) 2?rr \J dBr (x' ) n I 8vrr 



x The condition R > 1 could be replaced by R > Ro for any Ro > at the expense of a constant c depending 
on Rq. 
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Integrating over T, using \T\ > h(r — 1)R, we have 



dr 



dr > r R 
T T J T R-i(T-l)R r 



log 1 



and thus 



r- 1 
2r 



D(x ,R) 



[ \E Un \ 2 > cD(x ,R) 2 mm( 1 

Jb tR {x' )\b p V ^||m ||L°o 



for some c > depending only on | mo| l°° hence on V. Inserting into (4.1), we are led to 



f 9 Vn > -C(||mo|| £ cc + 1)U! + cD(x , R? min (l, ^ X °\ R) 
Jb 2R {x' ) V tt\\mo\\L°° 

The case D(xq,R) < is essentially analogous. 



□ 



We now proceed to the proof of Theorem pi starting with (1.36). If R > Rq and 



D(x' ,R)\ > rjR then from Proposition 6.5 and using the fact — from Proposition 



4.1 



that g Un is positive outside Uf =1 B(x' i ,C) and that g Un > —C everywhere, we deduce from 



(6.3) and (1.22), (1.23) that 



(6.11) 



fflK) > \ {-CR + cmmtf^^R 2 ) +2 J Qdv n . 



Inserting into (3.4) we find 



\D(x' ,R)\ >rjR) < -^exp(CPR-cf3mm(7] 2 ,7] 3 )R 2 ) f e~ nP ^ dVn dxi . . . dx n . 

Kn J 



Then, using the lower bound (3.6) and Lemma 6.3 we deduce that if ft > /3o and n is large 
enough depending on /?o then 

logF^(\D(x' ,R)\ >t)R) <-c(3mm(ri 2 ,T] 3 )R 2 + CpR + Cn(3 + Cn, 



where c, C > depend only on V. Thus (1.36) is established 



We next turn to ( |1.38 ). Arguing as above, from (6.3) we have F n (y n ) > —C + 2 f Qdv n - 
Splitting 2 J C,dv n as J Qdv n + J C,dv n , inserting into (3.4) and using (3.6) we are led to 



where C depends only on V. Then, using Lemma 6.3 we deduce (1.38). 



We finish with (1.37). Inserting the result of Lemma 6.4 into (3.4), we have, if I is an 



interval of width R/n 

F^(\\u n -n^\\ w ^, q{I) >C q rjVn(l + R 2 /n 2 )^^ < j d Xl . . . dx n . 



Arguing as before and rearranging terms yields (1.37). 
This concludes the proof of Theorem [6j 
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